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Abstract 

Orthogonal frequency division multiplexing (OFDM) has been considered for visible light communication 
(VLC); thanks to its ability to boost data rates as well as its robustness against frequency-selective fading channels. A 
major disadvantage of OFDM is the large dynamic range of its time-domain waveforms, making OFDM vulnerable 
to nonlinearity of light emitting diodes. DC-biased optical OFDM (DCO-OFDM) and asymmetrically clipped 
optical OFDM (ACO-OFDM) are two popular OFDM techniques developed for the VLC. In this article, we will 
analyze the performance of the DCO-OFDM and ACO-OFDM signals in terms of error vector magnitude (EVM), 
signal-to-distortion ratio (SDR), and achievable data rates under both average optical power and dynamic optical 
power constraints. EVM is a commonly used metric to characterize distortions. We will describe an approach to 
numerically calculate the EVM for DCO-OFDM and ACO-OFDM. We will derive the optimum biasing ratio in 
the sense of minimizing EVM for DCO-OFDM. In addition, we will formulate the EVM minimization problem as 
a convex linear optimization problem and obtain an EVM lower bound against which to compare the DCO-OFDM 
and ACO-OFDM techniques. We will prove that the ACO-OFDM can achieve the lower bound. Average optical 
power and dynamic optical power are two main constraints in VLC. We will derive the achievable data rates under 
these two constraints for both additive white Gaussian noise channel and frequency-selective channel. We will 
compare the performance of DCO-OFDM and ACO-OFDM under different power constraint scenarios. 

Index Terms 

Orthogonal frequency division multiplexing (OFDM), Visible light communications (VLC), DC-biased optical 
OFDM (DCO-OFDM), Asymmetrically clipped optical OFDM (ACO-OFDM), Error vector magnitude (EVM), 
Achievable data rate, Clipping 



I. Introduction 

With rapidly growing wireless data demand and the saturation of radio frequency (RF) spectrum, visible 
light communication (VLC) [|2) _ @ nas become a promising candidate to complement conventional RF 
communication, especially for indoor and medium range data transmission. VLC uses white light emitting 
diodes (LEDs) which already provide illumination and are quickly becoming the dominant lighting source 
to transmit data. At the receiving end, a photo diode (PD) or an image sensor is used as light detector. VLC 
has many advantages including low-cost front-ends, energy-efficient transmission, huge (THz) bandwidth, 
no electromagnetic interference, no eye safety constraints like infrared, etc. [6]. In VLC, simple and low- 
cost intensity modulation and direct detection (IM/DD) techniques are employed, which means that only 
the signal intensity is modulated and there is no phase information. At the transmitter, the white LED 
converts the amplitude of the electrical signal to the intensity of the optical signal, while at the receiver, 
the PD or image sensor generates the electrical signal proportional to the intensity of the received optical 
signal. The IM/DD requires that the electric signal must be real-valued and unipolar (positive-valued). 

Recently, orthogonal frequency division multiplexing (OFDM) has been considered for VLC; thanks 
to its ability to boost data rates and efficiently combat inter-symbol interference ||6)-[11|. To ensure that 
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the OFDM time-domina signal is real-valued, Hermitian symmetry condition must be satisfied in the 
frequency-domain. Three methods have been discussed in the literature for creating real-valued unipolar 
OFDM signal for VLC. 

(1) DC-biased optical OFDM (DCO-OFDM)— adding a DC bias to the original signal g, (5). 



(2) Asymmetrically clipped optical OFDM (ACO-OFDM) — only mapping the data to the odd subcarriers 
and clipping the negative parts without information loss |9); 



(3) Flip-OFDM — transmitting positive and negative parts in two consecutive unipolar symbols [11|. 



One disadvantage of OFDM is its high peak-to-average-power ratio (PAPR) due to the summation over 



a large number of terms [13|. The high PAPR or dynamic range of OFDM makes it very sensitive to 



nonlinear distortions. In VLC, the LED is the main source of nonlinearity. The nonlinear characteristics 



of LED can be compensated by digital pre-distortion (DPD) [14|, but the dynamic range of any physical 
device is still limited. The input signal outside this range will be clipped. A number of papers p"5|-[[T8| 
have studied the clipping effects on the RF OFDM signals. However, clipping in the VLC system has 
two important differences: (i) the RF baseband signal is complex- valued whereas time-domain signals in 
the VLC system are real-valued; (ii) the main power limitation for VLC is average optical power and 
dynamic optical power, rather than average electrical power and peak power as in RF communication. 
Therefore, most of the theory and analyses developed for RF OFDM are not directly applicable to optical 



OFDM. A number of papers [[14J, [[T9j-[(2T] have analyzed the LED nonlinearity on DCO-OFDM and 



ACO-OFDM and compared their bit error rate, power efficiency, bandwidth efficiency, etc. 

In this article, we will investigate the performance of DCO-OFDM and ACO-OFDM signals in terms 
of error vector magnitude (EVM), signal-to-distortion ratio (SDR), and achievable data rates. EVM is 



a frequently used performance metric in modern communication standards. In [14|, [20|, the EVM is 
measured by simulations for varying power back-off and biasing levels. In this article, we will describe 
an approach to numerically calculate the EVM for DCO-OFDM and ACO-OFDM, and derive the op- 
timum biasing ratio for DCO-OFDM. We will formulate the EVM minimization problem as a convex 
linear optimization problem and obtain an EVM lower bound. In contrast to p2| which investigated the 
achievable data rates for ACO-OFDM with only average optical power limitation, we will derive the 
achievable data rates subject to both the average optical power and dynamic optical power constraints. 
We will first derive the SDR for a given data-bearing subcarrier based on the Bussgang's theory. Upon 
the SDR analysis, we will derive the achievable data rates for additive white Gaussian noise (AWGN) 
channel and frequency-selective channel. Finally, we will compare the performance of two optical OFDM 
techniques. 

II. System model 

The system model discussed in this study is depicted in Figure 1 . In an OFDM system, a discrete time- 
domain signal x = [x[0], x[l], . . . , x[N — 1]] is generated by applying the inverse DFT (IDFT) operation 
to a frequency-domain signal X = [Xq, X\, . . . , X N _i] as 



x[n\ = IDFT(X fc ) = -= V X k exp(j2nkn/N), < n < N - 1, (1) 
v N k =0 

where j = and N are the size of IDFT, assumed to be an even number in this article. In a 

VLC system using LED, the IM/DD schemes require that the electric signal be real-valued and unipolar 
(positive- valued). According to the property of IDFT, a real- valued time-domain signal x[n] corresponds 
to a frequency-domain signal X k that is Hermitian symmetric, i.e., 

X k = X* N _ k , l<k<N-l, (2) 
X k e R, k = 0,N/2, 



where * denotes complex conjugate. 
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Fig. 1. OFDM system model in visible light communication 



In Figure 1, y[n] is obtained from x[n] after both a clipping and a biasing operation are implemented. 
The resulting signal, y[n], is non- negative (i.e., y[n) > 0) and has a limited dynamic range. In a VLC 
system, the light emitted is used for illumination and communication simultaneously. The intensity of 
light emitted by the LED is proportional to y[n), while the electrical power is proportional to y 2 [n]. In 
the optical communication literature, the average optical power of the LED input signal y[n) is defined as 

O y ± S{y[n}}, (3) 

where £[■] denotes statistical expectation. Usually, the VLC system operates under some average optical 
power constraint P A , i.e., 

O y < P A . (4) 

This constraint is in place for two reasons: (i) the system power consumption needs to be kept under a 
certain limit, (ii) the system should still be able to communicate even under dim illumination conditions. 
The VLC system is further limited by the dynamic range of the LED. In this article, we assume that the 
DPD has perfectly linearized the LED between the interval [Pl,Ph], where Pl is the turn-on voltage 
(TOV) for the LED. If the TOV is provided by an analog module at the LED and we assume that LED 
is already turned on, the linear range for the input signal is [0, Ph — Pl}- We define the dynamic optical 
power of y[n] as 

G y = max (y[n)) — min (y[n]) . (5) 
G y should be constrained by P H — Pl as 

G y < P H - Pl- (6) 

Moreover, y[n] must be non-negative, i.e., y[n] > 0. 

According to the Central Limit Theorem, x[n] is approximately Gaussian distributed with zero mean 
and variance cr 2 with probability density function (pdf): 

p(x[n] = z) = -S (-) , (7) 



cr Vcr, 

where Six) = -7=e~5 x ' 2 is the pdf of the standard Gaussian distribution. As a result, the time-domain 
OFDM signal x\n\ tends to occupy a large dynamic range and is bipolar. In order to fit into the dynamic 
range of the LED, clipping is often necessary, i.e., 

c u , x [tz] ^> c u 

x[n] = ^ x[n], Q < x[n] < c u (8) 
c h x[n] < q 
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where c u denotes the upper clipping level, and c% denotes the lower clipping level. In order for the LED 
input y[n) to be non-negative, we may need to add a DC bias B to the clipped signal x[n] to obtain 

y [n\ = x[n\ +B, 0<n<N-l. (9) 

For y[n] > 0, we need B = — cj. 

To facilitate the analysis, we define the clipping ratio 7 and the biasing ratio ^ as 

A (c u - Cl )/2 
a 

a B —a 
^ = ^^ = ^- (11) 

Thus, the upper and lower clipping levels can be written as 

C = 2<7 7 (1-?), (12) 

q = -2<x 7 g. (13) 

The ratios 7 and q can be adjusted independently causing c u and q to vary. 

Clipping in the time-domain gives rise to distortions on all subcarriers in the frequency domain. On the 
other hand, DC-bias only affects the DC component in the frequency-domain. The clipped and DC-biased 
signal y[n] is then converted into analog signal and subsequently modulate the intensity of the LED. At 
the receiver, the photodiode, or the image sensor, converts the received optical signal to electrical signal 
and transforms it to digital form. The received sample can be expressed as 

r[n] = (x[n] + B) <g) h[n] + w[n], (14) 

where h[n] is the impulse response of the wireless optical channel, w[n] is AWGN, and ® denotes 
convolution. By taking the DFT of Equation (14), we can obtain the received data on the A;th subcarrier 
as 

R k = X k H k + W k , k + 0, (15) 

where H k is the channel frequency response on the A;th subcarrier. 

Based on the subcarrier arrangement, DC-biasing, or transmission scheme, several optical OFDM 
techniques have been proposed in the literature. In this article, we will focus on the performance analysis 
of two widely studied optical OFDM techniques, namely, DCO-OFDM and ACO-OFDM. In the following, 
we shall use superscripts ( D ) and ( A > to indicate DCO-OFDM and ACO-OFDM, respectively. 

In DCO-OFDM, subcarriers of the frequency-domain signal X^ D ) are arranged as 

X ( D ) = [0 X[ D) X { 2 D) ... X^l, ... X< D) Xl {D) ] (16) 



where the 0th and AT/ 2th subcarriers are null (do not carry data). Equation (16) reveals Hermitian symmetry 
with respect to k = N/2. Let /Q denote the set of data-carrying subcarriers with cardinality The set 
of data-carrying subcarriers for DCO-OFDM is icf ] = {1,2,..., N/2 - 1, N/2 + 1, . . . , N - 2, N - 1} 
and \IO d D '\ = N — 2. The time-domain signal x^ D ^[n] can be obtained as 

2 N,2 ~ 1 / x 
x iD) [n] = -= (JK( X 1 D) ) cos(27rkn/N) - 3(X^ D) ) sm(2irkn/N)) , (17) 

fc=i 

which is real-valued. In DCO-OFDM, we first obtain a clipped signal x^ [n] similar to the procedure in 
([8]), and then add DC-bias B = — q to obtain the LED input signal 

x (D) [n] > c u 

y^[n]=x (D) [n}+B= { x {D ^[n}-c h c x < x^[n] < c u (18) 
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In the frequency domain, 



Yt D) =X { k D) + C k , Wk^O, (19) 



where C k is clipping noise on the kth subcarrier. 

In ACO-OFDM, only odd subcarriers of the frequency-domain signal X^ A ) carry data 

X( A ) = [0 X[ A) xi A) ... X {A ] 2 _ X X*$_ x ... X; [A) X* (A) ],(20) 



and X^ A ) meets the Hermitian symmetry condition (J2J). The set of data-carrying subcarriers for ACO- 
OFDM is K. (A) = {1,3, ...,JV - 1} and \)C d A) \ = N/2. Thus, the time-domain signal x^[n] can be 



obtained as 

JV/4-1 



(27r(2g + l)n/N) - 3(X^) sin (27r(2g + l)n/N)) , (21) 

g=0 

which is real-valued. It follows easily that x^ [n] satisfies the following negative half symmetry condition: 

x {A) [n + N/2] = -x (A) [n], n = 0, 1, . . . , N/2 - 1. (22) 

Denote by z[n] a generic discrete-time signal that satisfies z[n + N/2] = —z[n], n = 0, 1, . . . ,N/2 — 1 
and by z[n\ its clipped version where the negative values are removed, i.e., 

m = {M, z[n]>0, 23 
11 1 0, otherwise. 



It was proved in p3| that in the frequency-domain 

1 



2 Z k , Vfcodd. (24) 



In ACO-OFDM, we obtain the LED input signal y^[n] via 



y ( A )[ n ] = <j x W[ n ], < x^[n] < c u (25) 
0, x^[n]<0 



Equation (25 ) can be regarded as a 2-step clipping process, whereby we first remove those negative values 
in a^[n], and then replace those a;^ A )[n] values that exceed c u by c u . Since a^[n] satisfies (22), we 
infer based on ( [24] ) that 

yt A) = \x[ A) + C k , Vfc, (26) 

where C k is clipping noise on the A;th subcarrier in the frequency-domain. For ACO-OFDM, no DC-biasing 
is necessary and thus the biasing ratio q — 0. 

As an example, suppose that we need to transmit a sequence of eight quadrature phase-shift keying 
(QPSK) symbols. Table 1 shows the subcarrier arrangement for DCO-OFDM, whereas Table 2 shows the 
subcarrier arrangement for ACO-OFDM. The time-main signals x^ D ^[n] and x^[n] and the corresponding 
LED input signals y^[n] and are shown in 

Figure 2. We see that in x^[n], the last 16 values are a repetition of the first 16 values but with 
the opposite sign. It takes ACO-OFDM more bandwidth than DCO-OFDM to transmit the same message, 
although ACO-OFDM is less demanding in terms of dynamic range requirement of the LED and power 
consumption. 
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Fig. 2. An example of x <D) [n], y {D) [n], x {A) [n] and y {A) [n] to convey a sequence 8 QPSK symbols. For DCO-OFDM, 7 
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III. EVM ANALYSIS 

EVM is a figure-of-merit for distortions. Let X f = [Xl,X\, . . . , Jf]^] denote the iV-length DFT of 
the modified time-domain signal x^. EVM can be defined as 



C(x«,xt) 



£ 




2 


£ 




2 





(27) 



where XW = [X^\x\ r) 



(r) 

, -Xjy-i-J denotes the reference constellation. For DCO-OFDM, X k 
{D) . For ACO-OFDM, X^ = \x[ A) for k e 4 A) • 



for k e K, 



(r) 
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A. EVM calculation 

In DCO-OFDM, clipping in the time-domain generates distortions on all the subcarriers. We denote 
the clipping error power by = ^2 kelc (D) £[\X^ D ' — X^l 2 ]. Since the sum distortion power on the 
Oth and iV/2th subcarriers is small relative to the total distortion power of iV subcarriers, according to 
the Parseval's theorem, we can approximate Pffi as 

= £ 8 l\ x l D) ~ x { k D) \ 2 } (28) 

Af-1 

« ^S[\x iD) [n] -x {D) [n}\ 2 } 

n=0 

= N (J°°(z - c u fU (*) dz + jH (z - Cl fU Q dz^ 

= Na 2 (l + 4 7 2 (1 - q) 2 - 2 7 (1 - 00(27(1 " 0) - 27^(27?) 

-$(2 7 (1 - 0) - 4 7 2 (1 - 2< ^(27(1 - 0) + H ~ 270 + 4tV$( - 27OJ , 
where = f*^ <j>(t)dt. Thus, we obtain the EVM for the DCO-OFDM scheme as 



1 + 4 7 2 (1 - o 2 - 27(1 - 00(27(1 - 0) - 2 7 ^(2 7 

\ 1/2 

-$(2 T (l-0) -4 7 2 (l-0 2 $(2 7 (l-0) +$(-2 7 + 4 7 \ 2 $(-270 j • 

To find the optimum biasing ratio q*, we take the first-order partial derivative and the second-order partial 
derivative of Pffi with respect to the biasing ratio q 



Na 2 ^4 7 0(2 7 (1 - 0) - 470(2 T ?) (30) 
-8 7 2 (1 - 0$(27fc - 1)) + 8 7 2 ^( - 27OJ > 

iVo- 2 ^ 2 ^^?- 1)) +8 7 2 $( -2 7 0^- (31) 

We can see that if q = 0.5, dP^'/dq = 0. The second-order partial derivative d 2 P^P /dq 2 > for all 
q. Hence, if q < 0.5, dP^ /dq < 0. If q > 0.5, dPffi/dq > 0. Therefore, q* = 0.5 is the optimum 
biasing ratio which minimizes Py# . By substituting q* into Equation (29) we obtain the EVM for the 
DCO-OFDM scheme at the optimum biasing ratio as 



dq 



d 2p(D) 



dq 2 



' p(D) 

€S3 = 1/7^ < 32 > 



= v/2(l + - 27^(7). 

Remark 1: (i) <;* = 0.5 is the optimum biasing ratio for DCO-OFDM, regardless of the clipping ratio, 
(ii) When q = 0.5, we infer that c u = — q, i.e., when the x^[n] waveform is symmetrically clipped at 
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the negative and positive tails, the clipping error power is always less than that when the two tails are 
asymmetrically clipped (i.e., when c u ^ ci or when q ^ 0.5). 

Denote by e[n], n = 0, 1, . . . , iV — 1 a generic discrete-time signal with DFT E k , k = 0, 1, . . . , N — 1. 



When k is odd, E k can be written as 

„ N-l 



Ei 



kn 



^^e[n]exp (-j2tt- 

v n=0 

£ e N ex P (-i27r^ J + ^ I 

n=0 /v n=l 



(33) 



1 



22 e[n + N/2] exp ( -J 2 7r— - j'/ett 
o ^ 



1 ^ Z' /bTi 

E ( e N - e[n + iV/2]) exp ( -j2n — 

* n=0 



7fc=U 

Let = 2g + 1, g = 0, 1, . . . , N/2 — 1, Equation (33 ) can be further written as 



E 



2q+l 



N/2-1 

E 

n=0 
N/2-1 



[en 



[ n + JV/2])exp^-j27r-^ 



•j2tt 



iV 



(34) 



1 I / 

Therefore, {E k = E 2q+ i\ql^ ) l can be viewed as the DFT coefficients of a new discrete-time sequence 

{exp (-j2n§ ) (e[n] - e[n + iV/2])}^ 1 . Applying the Parseval's theorem to {exp (-j27rf ) {e[n\ - e[n + iV/2])}; 
we obtain, 

iV/2-1 1 AT/2-1 



£ |£ fc | 2 = £ \E 2q+1 \* = \ (e[n]-e[n + N/2}) 2 . 

A;=odd q=0 n=0 



In ACO-OFDM, we denote the clipping error power by = V {A ) 8 



cording to (35), we can calculate P 7 as 



X [ k A) /2-Xi 



p(A) 

7 



E « 



(A) 



(A) 



(35) 



Ac- 



(36) 



N/2-1 

E* 

n=0 



x (A) [n] 



N/2] 



x (A) [n + N/2] 



N r 

j£ [(x {A) [n] - x {A) [n] + x (A) [n + iV/2])' 



N 

~4 J_ 



\Z + C u ) 



- I dz+ — 



(z - cv 



= yd 2 ( - 2 7 (2 7 ) + $ (-2 7 ) + 4 7 2 $ (-2 7 ) ) 
Then we obtain the EVM for the ACO-OFDM scheme as 



dz 



e 



(A) 

7 



p(A) 



'4P 



7,0 



iV(x 2 



(37) 



\ £ fceX; (A) £ /2 
v/-4 7 (2 7 ) + 2$ (-2 7 ) + 8 7 2 $ (-2 7 ). 
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B. Lower bound on the EVM 
Let us consider the setting 

x[n] =x[n] +c[n), 0<n<N-l, (38) 

where x[n] is the original signal, c[n] is a distortion signal, and the resulting x[n] is expected to have a 
limited dynamic range 

max(i[n]) — min(x[n]) < 2^/cr. (39) 



In ( |38| ), all quantities involved are real-valued. 

Clipping can produce one such x[n) signal, but there are other less straightforward algorithms that can 



generate other x[n] waveforms that also satisfy (39). 
In the frequency-domain, 

X k = X k + C k . (40) 

Since x[n], c[n], and x[n\ are all real- valued, X k , C k , and X k all should satisfy the Hermitian symmetry 
condition §2§. Therefore, c[n] has the form 



N/2-1 

c\n 

w'Y 

k=l 

We are interested in knowing the lowest possible EVM 



2 iV/2_1 / \ 1 1 

-= V MC k ) cos(27ckn/N) - 3(C fc ) sin(2vrA;n/A^) + -^C + -^C N/2 cosfvrn) (41) 
VN <H \ ) VN VN 





c k 


2 ] 




x k 


2 ] 



^=1 h " ' ' (42) 



among all such x[n] waveforms. Afterwards, we can compare the EVM from the DCO-OFDM and ACO- 
OFDM methods to get a sense of how far these algorithms are from being optimum (in the EVM sense). 
We formulate the following linear optimization problem: 



minimize \C k 

k&K d 

subject to max (x[n]) — min < 27a 

N/2-1 

x\n\ = x\n 



2 



2 / \ 

+ — ( K(C fc ) cos(2nkn/N) - Qf(C fc ) sin(27rA;n/V) J (43) 
k=l V / 



+ ^=C + ^=C N / 2 cos(7m), 0<n<N-l 



1 r + 1 

y/N VN 
Cq, Cn/2 G R 

When the distortion of each OFDM symbol is minimized by the above convex optimization approach, 
the corresponding EVM of x[n) (which is proportional to £ [J2k&K, d l^l 2 ]) serves as me lower bound 
for the given dynamic range 2ja. 

C. Optimality for ACO-OFDM 

In this section, we will prove that the ACO-OFDM scheme achieves the minimum EVM and thus is 
optimal in the EVM sense. 
For ACO-OFDM, let us write 

x [n] = x {A) [n] + c[n], (44) 
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where c[n] is the clipping noise to ensure that x[n] has a limited dynamic range as described in (39). In 
the frequency-domain, we have 

X k = X { k A) + C k (45) 



where the Xu^ subcarriers are laid out as in (20). According to (35 ), when k is odd, the objective function 



in (43) can be written as 



N/2-1 



fc=odd 



n=0 



The dynamic range constraints in problem (43) can be viewed as two constraints put together. 

max (x {A) [n] + c[n}) - min (x {A) [n] + c[n\) < 2^a, < n < N/2 - 1, 
max {x {A) [n] + c[n]) - min (x (A) [n] + c[n}) < 2^a, N/2 <n < N — 1. 
Since = —x^ A '\n — N/2] when N/2 < n < N — 1, Equation (48) can be further written as 



max (rc^ln] — c[n + N/2]) - min (x {A > [n] - c[n + N/2]) < 2^a, 0<n< N/2 - 1. (49) 



M) 



From Equations (46), (47), and (49), the problem (43) can be recast as 

N/2-1 



minimize 



E 

n=0 



c 2 fnl 



subject to max [x {A) [n] + c[n]) - min [x {A) [n] + c[n]) < 2^a, < n < iV/2 - 1 
which is equivalent to 

N/2-1 



minimize 



E 

n=0 



c 2 \n] 



subject to x (A) [n] + c[n] < 27a, < n < N/2 - 1 
x {A) [n] + c[n] > 0, < n < N/2 - 1 



In Appendix, we prove that the solution c*[n] to (51 ) yields 

( 27a, x^[n]> 2 7 a 
x( A )[n] = x {A) [n] + c*[n] = I x (A ^[n], < x^[n] < 2 1 a 

{ 0, J x( A )[n] < 

and thus the ACO-OFDM scheme is optimum in the EVM sense. 



(46) 

(47) 
(48) 



(50) 



(51) 



(52) 



IV. SDR ANALYSIS 

Based on the Bussgang's theorem p4| , any nonlinear function of x[n] can be decomposed into a scaled 
version of x[n] plus a distortion term d[n] that is uncorrelated with x[n]. For example, we can write 



x[n] = a ■ x[n] + d[n], n = 0, . . . , N — 1. 



(53) 



Let R xx [m] = £{x[n]x[n+m]} denote the auto-correlation function of x[n], and let R xy [m] = £{x[n]y[n + 
m]} denote the cross-correlation function between x[n] and y[n] at lag m. For any given m, the correlation 
functions satisfy 

R xd [m] = 0, (54) 



Rsx[m] = aR xx [m]. 



(55) 
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Thus, the scaling factor a can be calculated as 



a 



Rxx [0] 
Rxx [0] 



(56) 



o- 



oo 



xx ■ p(x)dx 



Let /(•) denote the function linking the original signal to the clipped signal, it is shown in p5J that the 
output auto-correlation function R xx [m] is related to the input auto-correlation function R xx [m] via 



Rxxim] = ^2 



X) 65 



R xx [m] 



where the coefficients 



-IV 



rf<[exp(-|^)] 
m ^ dx. 



(57) 



(58) 



The input auto-correlation function R xx [m] can be obtained from taking IDFT of the input power spectrum 
density (PSD) 



R xx [m] = IDFT{P Xjfc } m , 



m 



0. 



N-h 



(59) 



where P 



X,k 



£\\X k 



is the expected value of the power on the feth subcarrier before clipping. Then 
it is straightforward to calculate the output PSD by taking the DFT of the auto-correlation of the output 
signal: 



X,k 



DFT{i2 5a [m]} fc , fc = 0,...,iV-l. 



Taking the DFT of Equation (53), the data at the kth subcarrier are expressed as 

X k = DFT{a -x[n]} k + DFT{d{n}} k 
= a- X k + D k , k e K d . 



(60) 
(61) 



Here, we assume that D k is Gaussian distributed, which is the common assumption when N is large [ 16 1. 
The SDR at the A;th subcarrier is given by 

S[\a-X k \ 2 ] a 2 P x , k a 2 P x , k 

k e ltd, (62) 



SDR fc 



e a 



P D . k P Xk - a 2 P x „ 



where P D , k = £[\D k 



P x k — a 2 P Xjk is the average power of the distortion on the kth subcarrier. 



According to Equation ( p6[ ), we can obtain the scaling factor a as a function of the clipping ratio 7 
and the biasing ratio <j: 



a 



1 



a- 



xx ■ p(x)dx 
1 . (z 



(63) 



z 2 -<f) 
a \a 



dz + 



QZ-<f) 



dz + 



G 



00 J 
C u Z-( 

a 



dz 



$ (2 7 (1 - q)) - $ (-2 7 <0 . 
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Note that in (63) we have used Equations (13]) and (12) for q and c u . According to Equation ( [58] ), we 
can obtain the coefficient hi as a function oftne clipping ratio 7 and the biasing ratio 



o-0(2 7 O - a<j) (27(1 - 0) ~ 2(J 7^ (-27?) + 2(77(1 - 0$ (~2 7 (1 - ?)) , £ = 

cr$(27(l-?))-o-$(-27?), £=1 
I CX P (-27 2 ^ 2 )^e (£ _ 2) (-2 7? ) - J= cxp (-2 7 2 (1 - 2 )#e ( ,_ 2) (2 7 (1 -<;)), l>\ 



(64) 



where He n (t) = (— l) £ exp 



is the probabilists' Hermite polynomials [26 1 



V. Achievable data rate 



In VLC, average optical power and dynamic optical power are two main constraints. Recall from 
Equation we can obtain the average optical power of y[n] as 

O y = £{y[n}} (65) 
= £{x[n]} + B 

= [ z—(p ( — ") dz + c u [ —<p ( — ") <iz + q / —0 f — ) cfc — q 
i Ci o- Vcr/ J Cu a \aJ J-ocC 

= a U{2j<-) - (27(1 - 0) " 2 7^ (-27?) + 27(1 " (-27(1 " 0) + 2 7^ 

Let a 2 w = £{w 2 [n}} denote the power of AWGN w[n), we define the optical signal-to-noise ratio (OSNR) 

O 

(66) 



OSNR 



Recall from Equation ([5]), we can obtain the dynamic optical power of y[n] as 

G y = max (y[n\) - min (y[n\) = c u - q = 20-7. 

We define the dynamic signal-to-noise ratio (DSNR) as 

G„ 



DSNR 



<7„ 



(67) 



(68) 



Let ?7osnr = Pa/&w denote the OSNR constraint and t/dsnr = (Ph — Pl)/<?w denote the DSNR constraint, 
we have 

(69) 
(70) 



O" < ??OSNR 



(7,: 



O y /a 



° < ^DSNR 



Gy/tr' 



The maximum o/a w value can be obtained as 



a 



mm 



^OSNR 



0(2 7 ^) - (2 7 (1 - q)) - 2 7 ^$ (-27?) + 2 7 (1 - q)$ (-27(1 - q)) + 2 7 <T 



^7dsnr \ 
27 / 



(71) 



by substituting (65 ) and (67 ) into the right-hand side of (69) and (70), respectively. The ratio //dsnr/^osnr 



(P H — P L )/P A is determined by specific system requirements. 
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AWGN channel 



For AWGN channel, recall from Equation (15), the received data on the kth subcarrier can be expressed 



as 



R k = X k + W k = aX k + D k + W k , ke K d . 
The signal-to-noise-and-distortion ratio (SNDR) for the kth subcarrier is given by 



SNDR fc 



a 2 £{\X k \ 2 } 



£{\D k \i} + £{\W k \*} 

D X,k 



a 2 P, 



Pd.k + cri 



(72) 



(73) 



SDR, + a. 



w a 2 Px,k 

In this article, we assume the power is equally distributed on all data-carrying subcarriers, 

Na 2 



P 



X,k 



then Equation (73) is reduced to 



SNDR fc 



1 



SDR, 1 + 2* • 



2 Nop 



By substituting Equation (71) into (75), we obtain the reciprocal of SNDR at the kth subcarrier: 

(SNDRfe)- 1 = (SDR,,)- 1 + S 

iVcr 



( 



■ max 



</>(2 7 <r) - <P (2 7 (1 - ?)) - 2 7 cr$ (-2 7 <r) + 2 7 (1 - (-2 7 (1 - <;)) + 2 7? 

2 

^OSNR 



4 7 2 



2 

^DSNR 



V 



(74) 



(75) 



(76) 



Therefore, the achievable data rate, as a function of clipping ratio 7, q, //osnr, and ^dsnr, is given by 

2 v bit s 

R (7, ?, ^OSNR, ??DSNr) = ^T7 5^ l0g2 ^ + SNDR *0 



subcarrier 



(77) 



Frequency-selective channel 

In the presence of frequency-selective channel, the received data on the A;th subcarrier obey the following 
in the frequency-domain: 

R k = H k X k + W k = H k (aX k + D k ) + W k , k e K d . (78) 



In this article, we consider the ceiling bounce channel model [27 1 given by 



h(t) = H(0) 



(t + a) 



-At), 



(79) 



where H(0) is the gain constant, a = 12-^/11/23.0 and u(t) is the unit step function. D denotes the rms 
delay. From Equation f78|), the SNDR is given by 



SNDR fc (# fc 





\Hk 


2 a 


2 £{\Xk\ 2 } 




\Hk\ 


2 £{\D k 


2 } + £{\W k 


2 } 



(80) 
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With the assumption of equal power distribution, we can obtain the 1 / SNDR as 

(SNDR,)- 1 = (SDR,)- + 

( ( N 2 

0(27?) - <j> (2 7 (1 - <;)) - (-27?) + 2 7 (1 - ?)$ (-2 7 (1 - ?)) + 2 7? 



2 

^OSNR 



V 

The achievable data rate, in the presence of frequency- selective channel, is given by 

j bits 

ft (7, ?, ^OSNR, ^DSNR, H) = — ^ log 2 (1 + SNDR fc (# fc )) 



(82) 



4 7 2 



' 2 

^DSNR 



subcarrier 



(83) 



VI. Numerical results 

In this section, we show EVM simulation results and achievable data rates of clipped optical OFDM 
signals under various average optical power and dynamic optical power constraints. 



A. EVM simulation 

The EVM analyses for DCO-OFDM and ACO-OFDM are validated through computer simulations. In 
the simulations, we chose the number of subcarriers N = 512, and QPSK modulation. One thousand 
OFDM symbols were generated based on which we calculated the EVM. In order to experimentally 
determine the optimum biasing ratio for DCO-OFDM, we used biasing ratios ranging from 0.3 to 0.7 in 
step size of 0.02, and clipping ratios ranging from 5 to 9dB in step size of 1 dB. Their simulated and 
theoretical EVM curves are plotted in Figure 3. As expected, the minimum EVM was achieved when the 
biasing ratio was 0.5, regardless of the clipping ratio. This agrees with the analysis in "EVM calculation" 
section. Next, we compared the EVM for DCO-OFDM with biasing ratio 0.5, EVM for ACO-OFDM with 
biasing ratio 0, and their respective lower bounds. To obtain the lower bounds, we used CVX, a package 



for specifying and solving convex programs [J28J, to solve Equation (43J). The resulting EVM curves for 
DCO-OFDM are plotted in Figure 4. The resulting EVM curves for ACO-OFDM are plotted in Figure 
5. We see that the EVM for ACO-OFDM achieves its lower bound, thus corroborating the discussion in 
"Optimality for ACO-OFDM" section. For DCO-OFDM, the gap above the lower bound increases with 
the clipping ratio (i.e., with increasing dynamic range of the LED). This implies that there exists another 
(more complicated) way of mapping x D [n] into a limited dynamic range signal x[n] that can yield a lower 
EVM. 



Achievable data rates performance 

We now show achievable data rates of clipped OFDM signals under various average optical power and 
dynamic optical power constraints. The number of subcarriers was N = 512. For the frequency-selective 
channel, we chose the rms delay spread D = 10 ns and sampling frequency 100 MHz. The normalized 
frequency response for each subcarrier is shown in Figure 6. 

As examples, we chose t/osnr = 20 dB, t^dsnr = 32 dB, and AWGN channel. Figures 7 and 8 show 
the achievable data rate as a function of the clipping ratio and the biasing ratio for DCO-OFDM and 
ACO-OFDM, respectively. We see that for given 7/osnr and ??dsnr values, a pair of optimum clipping 
ratio 7* and optimum biasing ratio q$ exist that maximize the achievable data rate. It is worthwhile to 
point out that the optimum biasing ratio ^ is different from (recall that minimizes the EVM). If 
the system is only subject to the dynamic power constraint, ^ should be equal to <j*. If the dominant 
constraint is the average power, ^ should be less than or equal to <j* because reducing the biasing ratio 
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Fig. 3. EVM as a function of biasing ratio for DCO-OFDM with clipping ratio = 5, 6, 9 dB 
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Fig. 4. EVM as a function of the clipping ratio 7 for DCO-OFDM along with the EVM lower bound for a given dynamic range limit 27a. 
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Fig. 5. EVM as a function of the clipping ratio 7 for ACO-OFDM along with the EVM lower bound for a given dynamic range limit 27a. 



can make the signal average power lower. We can obtain the optimum clipping ratio and biasing ratio for 
given ?7osnr, ^dsnr by 

(7W) = argmax ^| 7?OSNRi7?DSNR (84) 

Figure 9a shows the optimal clipping ratio as a function of ??osnr for DCO-OFDM. Figure 9b shows the 
optimum biasing ratio as a function of //osnr for DCO-OFDM. Similar plots are shown as Figure 10a,b 
for ACO-OFDM. In all cases, t/osnr varied from to 25 dB in step size of 1 dB, t/dsnr/^osnr = 18 dB, 
and the channel was AWGN. The main observation is, with a lower average optical power constraint, the 
clipping ratio and the biasing ratio can be increased to achieve higher data rates. Intuitively, when t^osnr 
is large, the channel noise has little effect and the nonlinear distortion dominates. 

Next, we chose the ratio ^dsnr/^osnr from 6dB, 12 dB, and no t^dsnr constraints. For each pair of 
Vosnr, Vdsnr, AWGN channel, or frequency- selective channel, we can calculate the optimum clipping 
ratio 7* and biasing ratio q$ according to Equation (84) and the corresponding achievable data rates. 
Figures 11, 12, and 13 show the achievable data rates with optimal clipping ratio and biasing ratio for 
the case ^dsnr/^osnr = 6dB, ^dsnr/^osnr = 12 dB, and no t^dsnr constraint, respectively. We observe 
that the performance of ACO-OFDM and DCO-OFDM depends on the specific optical power constraints 
scenario. In general, DCO-OFDM outperforms ACO-OFDM for all the cases. With the increase of the 
ratio ^dsnrA/osnr* the average optical power becomes the dominant constraint. The ACO-OFDM moves 
closer to the DCO-OFDM. 

As seen in Figure 13, when there is no DSNR constraint and the OSNR constraint is large, the DCO- 
OFDM curve closely matches the ACO-OFDM curve. This is in contrast to the performance curves 
in Figures 11 and 12. The reason for the curve coincidence in Figure 13 is twofold. First, we have 
already discussed that the performance difference between DCO-OFDM and ACO-OFDM is less when the 
OSNR constraint dominates, which is the case for Figure 13. Second, the suddenness of the convergence 
of the two curves can be explained by the fact that with only OSNR constraint, there will be more 
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Fig. 6. Normalized frequency response for each subcarrier (rms delay spread D = 10 ns, 100 MHz sampling rate). 
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Fig. 7. Achievable data rate as a function of the clipping ratio and the biasing ratio for DCO-OFDM with ijosnr = 20 dB, t]dsnr = 32 
dB, and AWGN channel. 
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Fig. 8. Achievable data rate as a function of the clipping ratio and the biasing ratio for ACO-OFDM with t/osnr = 20 dB, tjdsn r : = 32 
dB, and AWGN channel. 



flexibility in the signal optimization to adjust the clipping ratio and biasing ratio to achieve the best 
performance. That means that the achievable data rates in the middle-OSNR region (5-22 dB) are improved 
significantly compared with Figures 1 1 and 12. However, for high-OSNR region (greater than 22 dB), since 
the nonlinear distortion is negligible, the improvement becomes less pronounced compared to Figures 1 1 
and 12. Therefore, the transition from the middle-OSNR region to the high-OSNR region will become 
sharper with only an OSNR constraint. 

VII. Conclusions 

In this article, we analyzed the performance of the DCO-OFDM and ACO-OFDM systems in terms 
of EVM, SDR, and achievable data rates under both the average optical power and dynamic optical 
power constraints. We numerically calculated the EVM and compared with the corresponding lower 
bound. Both the theory and the simulation results showed that ACO-OFDM can achieve the EVM 
lower bound. We derived the achievable data rates for AWGN channel as well as frequency-selective 
channel scenarios. We investigated the trade-off between the optical power constraint and distortion. We 
analyzed the optimum clipping ratio and biasing ratio and compared the performance of two optical 
OFDM techniques. Numerical results showed that DCO-OFDM outperforms the ACO-OFDM for all the 
optical power constraint scenarios. 
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^OSNR ( dB ) 



Fig. 9. Optimal clipping ratio and biasing ratio of DCO-OFDM for tjosnr = 0, 1, 25 dB (in step size of 1 dB), t/dsnr/tiosnr = 18 
dB, and AWGN channel. 



Appendix 



Proof that c* [n] 



x 



(A) 



71 \ 



X 



(A) 



n] is optimum for Equation (51) 



Denote by w(c[n]) the objective function for the problem in (51 ): 

N/2-1 

u(c[n}) = ( C M) 2 , 



(85) 



n=0 



20 



(a) 




^OSNR < dB ) 



Fig. 10. Optimal clipping ratio and biasing ratio of ACO-OFDM for t/osnr ~ 0, 1, 25 dB (in step size of 1 dB), t/dsnr/tiosnr = 18 
dB, and AWGN channel. 



and denote by gi(c[n}) the ith constraint function for the problem in (51 ) 

9i(c[n}) = 



x (A ^[i] +c[i] -27a, 
-xW[i -N/2] -c[i-N/2], 

Let Hi denote the ith Kuhn-Tucker (KT) multiplier. We inter that 

Vtz(c[n]) = [2c[0],2c[l],... > 2c[iV/2-l]], 



< i < f - 1 
f <i <N-1 



N-l 



^ ^9i(c[n]) = [fl - fl N/2 , Hi - HN/2+1, VN/2-1 - Vn-i] ■ 



(86) 

(87) 
(88) 



t=0 
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Fig. 11. Achievable data rate with optimal clipping ratio and optimal biasing ratio for t/osnr = 0, 1, 25 dB (in step size of 1 dB), and 
Vdsnr/tiosnr = 6 dB. 



Next, we prove that c*[n] = x^[n] — x( A ">[n} satisfies the KT conditions [ |29| . 
Stationarity 



N-l 



Vu(c[n}) + J2 /^V</i(c[n]) = 0. 



j=0 



Primary feasibility 
Dual feasibility 
Complementary slackness 



Hi>0, Vz = 0, 1,...,N- 1. 
9i(c[n])<0, Vi = 0,l,...,iV-l. 

fi igi {c[n]) = 0, Vz = 0,l,...,jV-l. 



Substituting c*[n] into Equation (89), we obtain, 



Ik - lk+N/2 = -2x (A) [{\ + 2xW[i], i = 0, 1, . . . , — - 1. 



In order to satisfy all the other conditions <|90|>— <|92]>, we can choose m as follows 

(1) if x^[i\ = xW[i], 

fk — Hi+N/2 = 0; 

(2) if x^[i] > 2 7 a and x^[%] = 2-ya, 

fMi = 2x {A) [i], fi i+N/2 = 2x {A) [i]; 



(89) 

(90) 
(91) 

(92) 
(93) 

(94) 
(95) 
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Fig. 12. Achievable data rate with optimal clipping ratio and biasing ratio for ijosnr = 0,1,. ..,25 dB (in step size of 1 dB), and 
Vdsnr/tiosnr = 12 dB. 



(3) if x^[i] < and i< A )[i] = 0, 



Hi+N/2 



-2x^ A) [i 



X 



n 



(96) 

satisfy the KT 



Therefore, there exits constants //j (i = 0, 1, . . . , N — 1) that make c*[n] = x^ [n] 
conditions. It was shown in [ 30 1 that if the objective function and the constraint functions are continuously 
differentiable convex functions, KT conditions are sufficient for optimality. It is obvious that u and g are 
all continuously differentiable convex functions. Therefore, c*[n] is optimal for the minimization problem 
(|5T]>. 
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